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Abstract 

We perform an oscillation parameter-independent analysis of solar neu- 
trino flux measurements from which we predict the charged-current rate at 
SNO relative to Standard Solar Model to be R^§q = Q.?>btQ^^{'dAQtQ^^) 
for oscillations to active (sterile) neutrinos. By alternately considering the 

flux normalization fixed and free, we find that the flux measured by 
Super-Kamiokande (SK) not being a result of oscillations is strongly disfa- 
vored for oscillations to active neutrinos. SNO will determine the best-fit 
value of the ^B fiux normalization (5 (equal to the neutral-current rate), 
without recourse to neutral-current measurements, from the derived rela- 
tion (3 = R^Sj ~ 5.85i?sK — 4.85i?g^i^o- Using a simple parameterization 
of the fraction of high, intermediate, and low energy solar neutrinos starting 
above resonance, we reproduce the results of global analyses to good accu- 
racy; we find that the LMA solution with a normal mass hierarchy is clearly 
favored. With /3 free, our analysis for oscillations to active neutrinos gives 
rN^^ = P = 1.34 ± 0.34, which corresponds to R^^q = 0.28 ± 0.07. 



I. INTRODUCTION 



Solar neutrino experiments measure an energy-dependent flux suppression relative 
to the Standard Solar Model (SSM) 0. This situation has existed for thirty years and 
every new solar neutrino experiment has confirmed the fiux-deficit. The best motivated 
explanation of this solar neutrino puzzle is that neutrinos are massive and undergo oscilla- 
tions. The energy dependence of the flux suppression singles out very specific regions in the 
space of parameters that govern the frequency and amplitude of neutrino oscillations. The 
solar-neutrino flux deficit can be accounted for by oscillations of electron neutrinos to mu 
and/or tau neutrinos or to sterile neutrinos that do not interact weakly. In the two-neutrino 
oscillation framework, oscillations into sterile neutrinos are excluded at the 95% C. L. 
by a comparison of the day and night spectra at Super-Kamiokande (Super-K) and the 
results of a global flux analysis; hence the favored explanation is z/g —>■ v^/v-r oscillations. 
For such oscillations, there are three regions (LMA, SMA and LOW) in the mass-squared 
difference 5rn? and vacuum amplitude sin^ 26 parameters]^. These solutions involve effects 
from coherent z/g scattering from matter in the Earth and the Sun. Of these regions, the 
SMA region is disfavored at the 95% C. L. because of the observed fiat energy spectrum and 
imperceptible day/night effect at Super-K 0. For this reason we drop further consideration 
of the SMA region. A fourth region, the VAC or Just-So solution, which is largely indepen- 
dent of matter effects, is also excluded at the 95% C. L. by the same considerations that 
disfavor the SMA solution. Thus, we are left with the LMA and LOW regions, both of which 
have large mixing. With a large measure of certainty the KamLAND |[T^ experiment will 
exclude or confirm the LMA region as a solution to the solar neutrino anomaly |]Tli|l2| . The 



SNO experiment |13| will be crucial in validating the above conclusions from the Super-K 



experiment |]T4|,|T5| 



In this Letter we first perform a simple neutrino oscillation-independent analysis (with 
SSM fluxes) of the solar neutrino data using the total rates at the ^^Cl and ^^Ga ||-|] ex- 



periments and Super-K, following the procedure proposed in Ref. [|1^]. We make predictions 
for the charged- current (CC) rate at the SNO experiment. Allowing the flux normaliza- 
tion to be free, we derive a relation for the neutral-current (NC) rate at SNO in terms of the 
CC rate at SNO in a model-independent way. Our analysis is suited to the LMA and LOW 
solutions for which the oscillations in matter are mainly adiabatic. We find the relative flux 
suppression of the high, intermediate, and low energy solar neutrinos, compared to the SSM, 
and how these suppressions depend on the normalization of the solar '^B neutrino flux. We 
then apply our analysis to the LMA and LOW solutions and approximately reproduce the 
results of more comprehensive fits. Finally, we predict the CC and NC rates at SNO for the 
'^B flux normalization obtained by imposing adiabatic constraints. 



*From a flux-independent analysis the best-fit values of {5m? eV^,sin^2^) for the three solu- 
tions are (4.2 x 10^^0.66) (LMA), (5.2 x 10-^2.2 x IQ-^) (SMA) and (7.6 x 10-^,0.97) (LOW). 
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II. MODEL-INDEPENDENT ANALYSIS 



Following the procedure of Ref. |]T^, we divide the solar neutrino spectra into three 
parts: high energy (consisting of and hep neutrinos), intermediate energy (^Be, pep, ^^O, 
and ^^N), and low energy {pp). For each class of solar neutrino experiment the fractional 
contribution without oscillations to the expected rate from each part of the spectrum can be 
calculated in the SSM (see Table 1). If i? is the measured rate divided by the SSM prediction 
for a given experiment, then with oscillations, 

Rci = 0.764Pff + O.236P7 , (1) 
i?Ga = 0.096Pff + O.359P7 + 0.545Pl , (2) 

where Ph, Pi, and Pl are the average survival probabilities for the high, intermediate, and 
low energy solar neutrinos. In Eqs. (0) and (^ we assume that each probability remains 
the same from experiment to experiment, which is justified since the differential event rate 
without oscillations for each part of the spectrum has approximately the same shape for all 
experiments [p!^] . 

The SNO experiment detects neutrinos with energy above 5 MeV primarily via the 
reaction Ue + d ^ p + p + e~ ; the predicted CC rate is 

-^SNO = ■ (3) 

The R formula for Super-K is different from SNO, even though both experiments are sensitive 
to only the high energy neutrinos, because Super-K detects active neutrinos by elastic 
scattering, z/^ + e~ — z/^ + e~ (z/^ denotes any of the active flavors), and that has neutral- 
current contributions. Using the NC/CC cross section ratio of 0.171 (for xw = 0.225), we 
have 

RsK = 0.829Pff + 0.171 (active) , (4) 
= Ph (sterile) . (5) 

The solar neutrino data are summarized in Table |T[ Since there are three probability 
unknowns and three data points, there is an unique solution for the Pj. The Super-K rate 
depends only on the high energy neutrinos, so Ph is directly determined by -Rsk- Then since 
Rci depends only on Pj and Ph, the value of Pj may be determined from i?ci- Finally, Pl 
may be determined from Rca- The best-fit values to the data in Table |I| are 

TABLE I. Fractional contribution of the high, intermediate and low energy neutrinos to the 



CI, Ga and Super-K signals without oscillations. The last column gives the normalization 
uncertainty for each part of the spectrum. 







37C1 


7iGa 


Super-K 
SNO 


Normalization 
Uncertainty 


High 


'^B, hep 


0.764 


0.096 


1.000 


18.0% 


Intermediate 


^Be, pep, 1^0, 


0.236 


0.359 


0.000 


11.6% 


Low 


pp 


0.000 


0.545 


0.000 


1.0% 
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Experiment 


data/SSM 




0.337 ±0.030 




0.584 ±0.039 


Super — K 


0.459 ±0.017 



TABLE II. Solar neutrino data ||T|-|5| expressed as the ratio R = data/SSM, including the 
experimental uncertainties. The ^^Ga number combines the results of the GALLEX, SAGE, and 
GNO experiments. 



Ph = 0.347 , Pi = 0.303 , Pl = 0.811 (active) . (6) 
Ph = 0.459 , Pj = -0.058 , Pl = 1.029 (sterile) , (7) 

The sterile solution lies somewhat outside of, although close to, the physical region. 

To determine the allowed regions in probability space, we include the experimental un- 
certainties in the measured values and the theoretical uncertainties of the high, intermediate, 
and low energy neutrino fluxes in the SSM (see Table |). We use the following expression 
for x^'- 

/ p _ iDth\2 

^ ^ + alP]j5l + b^Pf6] + 4P151 ' 

Here i runs over the three types of experiments ('^^Cl, ^^Ga, and Super-K), Ri and 5Ri are 
the central values and uncertainties of data/SSM, 5j is the normalization uncertainty of that 
part of the solar spectrum (j = H, I, L), and the theoretical prediction for Ri is 

i?f = aiPn + hPi ± aPL + di . (9) 

The coefficients in Eq. (|^) are given in Eqs. ([l|), (0), (P, and (^, and the normalization 
uncertainties are given in Table |. 

The best-fit values and 95.4% C.L. {2a) allowed regions (Ax^ < 8.02), for each of the two- 
dimensional subspaces of the three-dimensional probability space are shown in Figs. |I] and |^ 
for oscillations to active and to sterile neutrinos, respectively. Note that the uncertainty in P/ 



spans the entire physically allowed region for active neutrino oscillations. The Borexino |17 



and KamLAND experiments will target ''Be solar neutrinos and should be able to narrow 
this range significantly. 

SNO observes only high energy neutrinos. Our analysis with the SSM fiuxes yields the 
simple prediction 

pg^Q = Ph = 0.35 (active) (10) 
= 0.46 (sterile) (11) 

with a possible range 

0.23 < Pg^o < 0.62 (active) (12) 
0.30 < R^^o < 0.71 (sterile) (13) 

at 95.4% C.L.. Because the ranges in Eqs. (|1^) and (|1^) are partially non-overlapping, a 
differentiation between the active and sterile scenarios from CC data alone could be possible. 
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III. UNCERTAINTY OF THE FLUX 
A. flux normalization flxed by Super K 

One approach to the uncertainty of the flux is to assume that the neutrinos are not 
suppressed at all by oscillations (or other particle physics mechanisms), and hence that the 
Super-K experiment is providing a direct measurement of the ^B flux. Then the implications 
of the ^''Cl and ^^Ga data on Pj and Pl can be examined. 

This scenario can be modeled by fixing Pfj = 0.459, the value of data/SSM measured by 
Super-K. Then the the remaining two types of experiments can be evaluated by using 

Eq. (H), with 6h set equal to the fractional uncertainty in the Super-K measurement, i.e., 
6h = 0.017/0.459 = 0.037, and summing over i = I,L. There are two degrees of freedom. 
Pi and Pl, and two data points. The unique solution, with = 0, is 

= -0.058, Pl = 1.029. (14) 

While this best-fit solution lies outside the physical region, a solution with Pj = and 
Pl = 1 is nearly as good, giving a of 0.2. At 95.4% CL. we find that P/ < 0.30 and 
Pl > 0.74. Hence for this scenario to be a good description of the data, the intermediate 
energy neutrinos are strongly suppressed, while the low energy neutrinos are not greatly 
suppressed. Only a vacuum solution with 5m? ~ 6 x 10~^^ eV^ and large mixing can give 
such probabilities |T^; however, the global analysis of Ref. 0] shows that this solution is 



acceptable only at 99% C. L. for active neutrino oscillations (90% C. L. for oscillations to 
sterile neutrinos). If none of the low or intermediate energy neutrinos are suppressed, the 
is 91.2. Thus the scenario in which the ^B flux is not suppressed by oscillations and is 
being directly measured by Super-K is highly disfavored for oscillations to active neutrinos. 



B. Varying the ^B flux normalization 

The uncertainty in the ^B flux is relatively large (18% at the la level), and it is possible 
that the SSM does not give a good estimate of it. Rather than just include the ^B flux 
uncertainty in the calculation of x^? we can allow the *B flux normalization to be a free 
parameter in the fit. We define f3 to be the ^B flux relative to the SSM. For oscillations to 
active neutrinos, the neutral current rate (NC/SSM) at SNO is 

<NO = P (active) . (15) 

for oscillations to active neutrinos is determined by using Eq. (^ with the changes 6h ^ 
(the uncertainty in the ^B flux is to be determined by the fit), Ph — > PPh in Eqs. (|l]) and 
(I), and PsK ^ /3{0.829Ph + 0.171) in Eq. (|). 

Since there are now four parameters and only three data points, there is no longer a 
unique solution, but a family of solutions parameterized by j3. The values of Ph, Pi, and 
Pl that exactly reproduce the data are plotted versus (3 in Fig. |^. The figure shows that 
for no value of P are the Pj all equal (although two of the three could be equal), so that 
there must be an energy-dependent suppression compared to the SSM. Values of (3 < 1 (i.e., 
initial ^B flux less than the SSM) imply a greater suppression of the intermediate energy 
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oscillation probability while for /3 > 1 the high energy probability is more suppressed. 
Figure |^ represents the possible exact solutions of the solar neutrino puzzle against which 
particular models can be compared. The only caveat is that our analysis accounts only 
for the average rates for each part of the neutrino spectrum; other measurements, such 
as energy dependence within one part of the spectrum or the day/night asymmetry, may 
provide further constraints. 

The family of exact solutions give a prediction for SNO that depends only on (3 and -Rsk, 

RsNO = PPh = 1.21i?sK - 0.21/5 (16) 
= 0.55 - 0.21/3 {Rsk = 0.459) . (17) 

Thus the SNO CC measurement will select a particular best-fit value of the '^B fiux normal- 
ization and the NC rate, 

^^No = P = 5.85i?sK - 4.85i?g?o (18) 
= 2.69 - 4.85i?g^o (^SK = 0.459) . (19) 

Performing a fiux- independent analysis by varying Ph, Pi, Pl and /3, the predicted CC 
rate at SNO has the upper bound 

^^No < 0.5 (*B fiux free) (20) 

at 95.4% C.L.. The large allowed range of R^^q represents the variation of P from 0.4 to 
2.8 (see the plot of (3 versus pPn in Fig. H). The best-fit line with = (since there 
are four parameters and three constraints), and 95.4% C.L. (Ax^ < 9.70) allowed regions 
are shown in Fig. ^ The plots are made versus pPn to make the uncertainties in R^^q 
transparent. Once adiabatic constraints are included for the MSW solutions with large 
mixing, the number of free parameters is three, and the best-fit line collapses to a best-fit 



point marked by the cross; we will discuss this in Section |1VB . 

For sterile neutrinos, Ph is replaced by pPn in Eqs. (|I|), (^, and (^. There is again a 
family of solutions for the probabilities, but unlike the active case, the values of Pj and Pl 
are fixed, and are the same as those in the sterile solution with the ^B normalization fixed 
at unity, Eq. (|^. Only the value of Ph depends on /?, with pPn = 0.459. Thus we see that 
in the sterile case the intermediate energy neutrinos must be strongly suppressed regardless 
of the value of the ^B normalization. This shows why the sterile case requires the SMA 
solar solution, which has a strong suppression of the intermediate energy neutrinos. Since 
-^SNO ~ = 0.459, the best fit solutions for sterile neutrinos give the same prediction for 
-^SNO f*-*^ value of the ®B fiux normalization. 

Because the SMA solution seems increasingly unlikely 0,0] for either active or sterile 
oscillations, we consider only active oscillations in the LMA and LOW solutions. 



IV. MSW SOLUTIONS WITH LARGE MIXING 

The consistency of particular neutrino oscillation solutions can be tested using the proba- 
bilities Ph, Pi, and Pl determined in the previous sections. In MSW solutions the oscillation 



probability is |19 
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P(l/e ^ I^e) = 2 + ( 2 - COS 20 , (21) 

where Om is the effective mixing angle in matter at the creation point of the electron neutrino 
and Pc is the Landau-Zener probability for crossing from the upper to the lower eigenstate 
as the neutrino propagates through the matter in the sun. We consider only the LMA and 
LOW solutions, whose suppression of neutrinos has very little energy dependence, in 
agreement with the measured Super-K spectrum. 



A. Adiabatic solutions with B flux from SSM 

In a typical MSW solution with large mixing, all of the solar neutrinos propagate adi- 
abatically, which implies Pc = in Eq. (^) |2^. For neutrinos created in a region of the 
sun above the critical density for a resonance to occur and that start far above resonance, 
cos20A/ = —1, which implies the oscillation probability is 

P(z/e z/e) = 1(1 - cos2^) = sin'^. (22) 

For neutrinos that start well below resonance, cos 26' a/ = cos 26* and 

P(i/e ^ I/,) = 1(1 + cos' 2^) = 1 - 1 sin' 26 . (23) 

For each part of the energy spectrum, we define fj {j = H, I, L) as the fraction of those 
neutrinos that are created above resonance; then 

P,- = /, sin' e + il- /,)(1 - 1 sin' 29) . (24) 

In Eq. ( |24D we have implicitly assumed that a negligible fraction of neutrinos are created 
near the resonance. Although this is not strictly true for a neutrino of a given energy, when 
averaged over the entire spectrum Eq. (^41) provides a good approximation to the overall 
average probability P, . 

The condition for which neutrinos are created above resonance is 

W cos 29 < 2V2GFNeE = (1.52 x 10^^ eV^)N^E , (25) 

where is the electron number density in the Sun in units of A^/cm^ and E is the neutrino 
energy in MeV; for a typical high or intermediate energy neutrino in the sun, A'g — 90 [§. 
Thus, for a given 5m' and 9, Eq. (25) defines the critical neutrino energy above which 



neutrinos are created above resonance. 

The resonance condition depends directly on the neutrino energy. For a given Jm' all 
of the neutrinos above a certain critical energy will be created above resonance. Since 
neutrinos in one part of the spectrum cannot be above resonance until all of the neutrinos 
with higher energy are also above resonance, we see immediately that MSW solutions must 
have /iT > // > /l- Furthermore, we can define three regimes: (i) some or all of the high 
energy neutrinos and none of the low or intermediate energy neutrinos are created above 
resonance (0 < < 1, // = /l = 0), (ii) all of the high energy neutrinos, some or all 
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of the intermediate energy neutrinos, and none of the low energy neutrinos are created 
above resonance {/h = 1, < // < 1, = 0), and (iii) all of the high and intermediate 
energy neutrinos and some or all of the low energy neutrinos are created above resonance 
{Ih = // = 1, < /l < 1). It should be noted that although there are a priori three fraction 
parameters {fn, fi, and fi) in an MSW solution, // cannot be nonzero unless /i^ = 1 and /l 
cannot be nonzero unless // = 1, so the fj are in fact equivalent to a one-parameter system 
with each value of Srri^ corresponding to an unique set of /-,-. By examining the expected 
differential rates versus neutrino energy for a given part of the spectrum, the values of the 
fj may be determined. For example, when sin^ 26' = 0.8 (cos 26' = 0.45), = 10~^ eV^ 
implies a critical energy of about 3 MeV, so that all of the high energy neutrinos and none 
of the low and intermediate energy neutrinos are created above resonance; this corresponds 
to /h = 1 and // = /l = 0. Thus, 9 and the fj are effectively two free parameters. 

In Fig. 1^ we show versus the fj for these three regimes for various values of sin^ 9, 
assuming the normalization is given by the SSM with the uncertainty in Table |. The 
best fit is 

sin^ 9 = 0.273 , fn = 0.874 , /, = = , (26) 

which in terms of probabilities, Eq. (p^ is 

Ph = 0.315 , Pi = 0.603 , Pl = 0.603 . (27) 

The x^/d-o.f. is 1.02/1, which corresponds to a goodness-of-fit of 31%. The associated 
oscillation amplitude is sin^ 29 = 0.79, which is almost exactly the value obtained in the 
flux-constrained global analysis of solar neutrino data in Ref. For this value of fn, 87% 
of the high energy neutrinos are created above resonance. An inspection of the ^^Cl and 
Super-K spectra shows that neutrinos with energies above about 7 MeV are created above 
resonance, which by Eq. (^) gives ~ 2 x 10~^ eV^, a value that corresponds to an LMA 
solution. However, from the fiux-dependent global analysis of Ref. |Q , this value of Srn^ is 
allowed only at the 99% C. L. which indicates that the day and night spectra and the full 
resonance treatment with unaveraged energy spectra are important in the global fit. 

For the LOW solution, in our approximation all of the neutrinos are created with 
cos29m = -1 in Eq. (|l|), so fu = fi = fi = I and Pjj = Pj = P^ = sm^9. This 
condition is far from the exact fit of Eq. (^). The best fit for Ph = Pi = Pl with the SSM 
^B fiux normalization (i.e., /5 = 1) is Pj = sin^ 6' = 0.517 with x^/d.o.f. = 10.6/2, which is 
excluded at the 99.5% C.L. Thus the LOW solution is disfavored. 

One might wonder why the values of Pj in Eq. (pTj) are not the same as those in Eq. (^ 
from the model-independent analysis. To understand this, we note that the situations in 
which the intermediate and low energy neutrinos can have different flux suppressions are 
cases (ii) and (iii) described above, both of which require //^ = 1. Working backwards 
from Eqs. (^ and (^) with /jy = 1 gives // = 1.22, which is unphysical. Thus there is 
either an inconsistency in the data from the different experiments or not all the SSM flux 
normalizations are correct. 
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B. Adiabatic solutions with free B flux normalization 

The MSW large angle solutions can also be tested when the normalization is allowed 
to vary, in which case the Pj are given by Fig. ^ as a function of the normalization factor 
13. The figure shows that for /3 < 1 (^B flux less than the SSM prediction) Ph and Pl 
are driven higher and Pi lower than the best fit for /3 = 1. This would be very hard to 
understand in the LMA solution, because the high energy neutrinos typically start above 
resonance (with P = sin^ 9) and the low and intermediate energy neutrinos typically start 
below resonance (with P = 1 — ^ sin^ 29); since sin^ 9 < 1 — ^ sin^ 29 for 9 < 7r/4, this implies 
Pfj < Pi < Pl- On the other hand, for (3 > 1, Pj and Pl are higher than P^, which is 
qualitatively consistent with the LMA solution. By imposing the condition Ph < Pi ^ Pl, 
required for an MSW solution with 9 < 7r/4, Fig. ^ shows that the LMA solution favors a 
'^B flux normalization in the range 1.04 < /3 < 1.48. 

We note that for an MSW solution with 9 > it/ 4, all of the neutrinos are created 
with cos 29m = — 1 in Eq. (|21|), so in our approximation fn = fi = II = ^ and Ph = 
Pi = Pl = sin^ 9. There is no value of P for which this occurs in the overall best fit (see 
Fig. hence this region of parameter space is disfavored, and a normal mass hierarchy 
{9 < 7r/4,(5m^ > 0) is selected. Similarly, as discussed above, all three probabilities are 
equal in the LOW solution in our approximation; the best fit for f3 and Ph = Pi = Pl free 
is /3 = 0.72 and Pj = sin^ 9 = 0.540 {9 > 7r/4), with xVd.o.f. = 9.16/1. Thus the LOW 
solution is approximately 3a from the overall best fit, similar to the findings of Ref. [p. 

In the above modelling with free ^B normalization, there are three free parameters for 
the large angle solutions: the vacuum mixing angle 9, the ^B normalization (3, and the 
fractions fj that determine which neutrinos start above resonance. For oscillations to active 
neutrinos, we find that the unique solution (with = 0) is 

sin^ 9 = 0.206 , (3 = 1.34 , = 1 , /, = 0.301 , /l = . (28) 

This solution is very close to the values sin^ = 0.21 and P = 1.31 found in the global 
analysis of Ref. , which also allowed the ^B flux be a free parameter. The fj of Eq. (p8[) 
imply that all of the high energy neutrinos and about 30% of the intermediate energy 
neutrinos are above resonance, which implies that the critical energy lies on the ^Be line at 
0.862 MeV; this translates into 6m^ = 2 x 10"^ eV^ for sin^ ^ = 0.206 (sin^ 20 = 0.654), 
close to the best-fit value of Ref. 0. Thus our simplified analysis reproduces the results of 
a comprehensive global fit to the data. There is no exact solution possible for oscillations to 
sterile neutrinos since that would require Pi < Ph, Pl, which is not allowed by the ordering 
Ph ^ Pi Pl', this explains why oscillations to sterile neutrinos are disfavored for MSW 
solutions with large mixing angles. 



The probabilities corresponding to Eq. (^81) are (using (|2J)) 

Ph = 0.206 , Pi = 0.532 , Pl = 0.673 . (29) 
The best-fit values of and uncertainties in /3 (R^m)) -^sno (^^e Fig. H) 

R^^Q = (3 = 1.34 lb 0.34 (30) 

and 
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R^^^ = (3Ph = 0.28 ± 0.07 . (31) 



The three values of Pj, and /5 from Eq. ( pSD represent a unique best-fit in Fig. ^ marked by 
a cross; the adiabatic constraints on the hne of best-fit solutions singles out one solution. 



V. SUMMARY 

By parameterizing the expectations for the three types of solar neutrino experiments 
(^^Cl, ^^Ga and u-e scattering), in terms of three average survival probabilities for the high, 
intermediate, and low energy solar neutrinos, we have determined a unique best fit assuming 
SSM fluxes. Accounting for the experimental and theoretical uncertainties, allowed regions 
in the probability space were found. Our analysis with the SSM fluxes yields a CC prediction 
for data/SSM for SNO of 

i?s% = 0.35^°:^^ (active) (32) 
= OAQtofi (sterile) (33) 

where the uncertainties are Icr. The prediction for oscillations to sterile neutrinos is fiux- 
independent. For some values of R^^q it could be possible to distinguish between the active 
and sterile scenarios without using with the neutral-current rate. 

A scenario in which the flux is not affected by oscillations and is assumed to be directly 
measured by Super-K is highly disfavored by the data for oscillations to active neutrinos. 
Allowing the normalization of the *B flux to be a free parameter, a family of solutions was 
found that depend on the flux normalization factor; ^B flux normalizations below the SSM 
imply that the intermediate energy neutrino contribution must be more suppressed, while 
normalizations above the SSM imply that the high energy neutrino contributions are more 
suppressed. The family of best-fit probabilities of Fig. ^ give a /3-dependent prediction for 
the best-fit CC rate at SNO: 



RsNO = PPh = 1.21i?sK - 0.21/5 (34) 
= 0.55 - 0.21/3 (RsK = 0.459) . (35) 

This equation can be inverted so that the central value of SNO CC measurement determines 
the central value of the NC rate, 

^^NO = P = 5.85i?sK - 4.85i?g^o (36) 
= 2.69 - 4.85i?g^o (^SK = 0.459) . (37) 



Thus, it may be possible for SNO to obtain the ^B flux normalization without recourse to 
neutral-current measurements. 

By imposing adiabatic constraints on our probability parameterization with the ^B flux 
free, we found a unique solution where all of the high energy and 30% of the intermediate 
energy neutrinos are created above resonance. This solution has a best-fit *B flux normal- 
ization (3 = 1.34, and 

W = 2 X 10"^ eV^ , sin^ 26 = 0.65 , (38) 
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in the LMA region. These values are close to those found from more comprehensive global 
fits to solar data |^. The LOW solution and an inverted mass hierarchy with 6 < n/A and 
Sm'^ < are disfavored. From our best-fit, we predict 

P = R^^Q = 1.34 ± 0.34 (39) 
R^^o = 0.28 ± 0.07 (active, best fit fiux). (40) 
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FIG. 1. Best-fit values and 95.4% C.L. allowed regions for Pj versus Ph, Pl versus Pr, and 
Pl versus Pi, assuming oscillations to active neutrinos and SSM fluxes (/3 = 1). 
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FIG. 2. Same as Fig. ffl for oscillations to sterile neutrinos. 
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FIG. 3. Probabilities Ph, Pi, and Pl for active neutrino oscillations that exactly reproduce the 
solar neutrino data, plotted versus /3, the ratio of the flux normalization to its SSM value. 
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FIG. 4. The 95.4% C.L. allowed regions from an analysis with free flux normalization /3, for 
Pj versus (3Ph, Pl versus fiPn, Pl versus Pi and j3 versus j3PHi assuming oscillations to active 
neutrinos. The line represents solutions with = 0. The cross marks the point selected once 
adiabatic constraints are imposed. 
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FIG. 5. for the LMA solution versus (i) fn, with fi = fL = 0, (ii) //, with fn = ^ and 
/l = 0, and (iii) fi, with fH = fi = ^, where fj is the fraction of j-type neutrinos {j = H, I, or 
L) created above resonance, for various values of sin^ 6. SSM fluxes are assumed (/3 = 1). 



16 




FIG. 6. The 95.4% C. L. and 68.3% C. L. allowed regions from an analysis with /?, and fj 
free. The best-fit point (iisNb' -^SNo) = (1-34,0.28), is marked with a cross. 
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